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We demonstrate that network models for wave mechanical systems with quenched disor- 
der cover the physics of mesoscopic electrons. The models are constructed as a network 
of random scattering matrices connecting incoming to outgoing wave amplitudes. The 
corresponding wave dynamics is given by a discrete unitary time evolution operator. We 
report on three different universality classes: two-dimensional, spinless, non-chiral elec- 
trons with (02NC) and without time reversal symmetry (U2NC), and two-dimensional, 
non-chiral electrons with time reversal symmetric spin-scattering (S2NC). We determine 
the phase diagram in the parameter space of scattering strengths. The 0/U2NC mod- 
els show strong localization. We find symmetry factors in localization lengths as well 
as multifractal exponents in agreement with theoretical predictions. The S2NC model 
displays a localization-delocalization transition. We determine the critical exponent of 
the localization length and the multifractal scaling exponent of the order parameter to 
be v ft: 2.4 and olq 2.18, respectively. 



In phase coherent disordered electron systems (mesoscopic systems) interference 
phenomena can lead to strong fluctuations in physical quantities and localization 
of wave functions. Localization-delocalization (LD) transitions tuned by system 
parameters can be described in terms of the scaling theoryEJ for broad distribution 
functions of physical quantities 13. In the present work we start from a modeling 
of mesoscopic systems by networks of unitary scattering matrices (see Fig [[]) . 
These models allow for a direct evaluation of physical quantities and of (quasi- 
energy) eigenvalues and eigenstates. Network models (NWMs) are paradigmatic for 
coherent waves in disordered media and do not refer to any particular dispersion 
relation. They are designed to cover essential symmetries and characteristic length 
scales, but are otherwise unspecific. For example, the wavelength can be identified 
with the lattice spacing which is, together with the wave velocity, set to unity. 
The elastic mean free path l e is determined by the average reflection properties of 
single scattering units. NWMs can rather easily be designed for special purines. 
For example, the situation of strong magnetic fields in 2D can be modeled QO by 
suppression of forward and backward scattering in each scattering unit of Fig. [I]. 




Figure 1: A two-dimensional network of scatterers (squares) connecting incoming and outgoing 
propagating wave amplitudes (bonds). 
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Figure 2: The parameter space of the U2NC network model: r 2 + t 2 < 1 where r and t denote 
transmission strengths of individual scatterers as shown in the inset. Vanishing deflection leads to 
decoupled ID chains, maximum transmission (reflection) corresponds to the metal (localization) 
fixed point, and maximum deflection corresponds to the quantum Hall fixed point. Systems 
with time reversal symmetry (02NC) have a restricted parameter space, r + t > 1, denoted as 
'orthogonal symmetry'. 

Consider a regular network of J\f sites and N bonds as shown in Fig. [|. The 
bonds carry propagating waves (incoming and outgoing) represented by complex 
amplitudes. On the sites unitary S'-matrices map incoming to outgoing amplitudes. 
The elements of each S'-matrix are (in general) random quantities taken from certain 
distribution functions, characterized by a few parameters, e.g. scattering strengths 
as shown in the inset of Fig. |^. Random phases (respecting symmetries) are attached 
to the amplitudes on the links. They simulate random distances between scatterers 
in realistic systems. The construction of a NWM is fixed by the choice of a type 
of S- matrix and a connectivity matrix C that describes how sites -am connected to 
each other. S and C define a unitary time evolution operator UuB that maps all 
incoming to outgoing bond amplitudes in one unit of time, 



Here the state tp is the vector of the N bond amplitudes. The eigenphases of U 
are appropriate objects for investigating local energy level statistics. The NWM 
depicted in Fig. |l| (with the scattering unit shown in the inset of Fig. [|) is designed 
to describe 2D disordered, non-interacting, spinless electrons in the absence of chiral 
fields. Its parameters are the transmission (reflection) strength t (r), and a deflection 
strength d (equal for left and right) obeying the constraint of unitarity, t 2 +r 2 +2cP = 
1. The parameter space is shown in Fig. g. The model is denoted as U2NC model. 
If the scattering is time reversal symmetric the parameter space is further restricted 
by r + t > 1 and the corresponding NWM is referred to as 02NC model. Under 
real-space renormalization the U2NC (02NC) model has three (two) fixed points: 
metal, localization, and quantum Hall fixed point □ (metal and localization fixed 
point). We have calculated localization lengths £(M) in quasi-lD strip geometries 
of width M, and the multifractal f(a) spectrum of eigenstates in a square geometry. 
It turns out that only the localization fixed point is attractive under renormalization 
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Figure 3: On the left normalized localization lengths A are shown as function of the inverse strip 
width 1 /M for different symmetries (/3 = 1,2) and different choices of the scattering strengths t 
(r was set to r = 1 — t) . On the right the multifractal spectra of states in a square geometry are 
shown for different choices of the scattering strengths and j3 = 2. The scaling exponent an is given 
as the maximum position of /(«). 



while the others are repulsive, i.e. all states will localize for large enough system 
sizes. The elastic scattering length is given as l e — (t 2 + d 2 )/(2r 2 + 2d 2 ) (for details 
seeB) and by the Einstein relation a diffusive classical conductance can be defined as 
o = ^e/(27r). As a quantitative test we compared £(Af) with the analytic expression 

£(M)=/ e (/?(2M) + 2-/3). (2) 

Here j3 = 1,2 is the usual symmetry index indicating presence (1) or absence (2) 
of time reversal symmetry. Equation (|J) is expected to be valid in strip geometries 
for which A = £(M)/M 3> 1. By finite size scaling this corresponds to a 2D 
metallic system in the weak localization regime. In fact, we find good agreement 
(with an uncertainty of 10%) already for moderate values of A > 1 (see Fig. ||) telling 
that the network model is able to describe disordered 2D electrons. As a second 
quantitative test we calculated the scaling exponenLao describing the scaling of the 
typical local density of states in a square geometryu. Such states correspond to the 

feak localization regime. In 2D ag can be related to A(M) by a conformal mapping 
as soon as A becomes independent of M (in practice this means At 3 for (3 = 2). 

A = ^Y)^ a °- 2 = Wio (3) 

This finding coincides with the analytic result0 obtained for large classical conduc- 
tance go and is consistent with our numerical results for an, the maximum positions 
in the multifractal /(a) spectra displayed in Fig. ||. The NWM describing 2D dis- 
ordered time reversal symmetric electrons in the presence of spin scattering (S2NC 
model) is based on the 02NC model til. In addition to the scattering parameters 
t,r of the 02NC model a new parameter, the spin scattering parameter s £ [0, 1], 
appears. It defines a spin scattering length ls(s) 6 [0, oof. The S2NC model shows 
true LD transitions in the parameter space (r, t) for non-zero values of s. The 
transitions can be identified by finite size scaling techniques based on the quantity 
A(M). The phase diagram obtained for maximum spin scattering strength s = 1 
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Figure 4: On the left the phase diagram of the S2NC model is shown for maximum spin-scattering 
strength s = 1. The grey area stands for the delocalized phase (D) separated from the localized 
phase (L). On the right a one-parameter scaling function is shown for the LD transition in the 
S2NC model. Data correspond to the logarithm of A as a function of ln(M/^ c (r)) where £ c (f) 
denotes a fitted correlation length as a function of the reflection strength r. 



is shown in Fig. |j. When crossing the phase boundary in the (p. t) plane, for a 
fixed value of s, the quantity A(M) follows one-parameter scaling^ as can be seen 
in Fig. ^. Analyzing the scaling function allowsJbr the determination of the critical 
exponent of the localization length, v w 2.4 til With the help of the conformal 
mapping relation (see Eq. (||)) the multifractal exponent can be obtained from the 
fixed point value of A and is olq « 2.18. 

In summary, we have shown that NWMs represent universality classes of disor- 
dered wave mechanical systems. In particular, the U/02NC models show Anderson 
localization. A quantitative analysis in the weak-localization regime is in reasonable 
agreement with known analytic results for disordered electron systems. The S2NC 
model describing spin-scattering exhibits localization-delocalization transitions and 
allows for a quantitative analysis of critical properties. 
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